The impact of non-LTE effects and granulation inhomogeneities on 
the derived iron and oxygen abundances in metal-poor halo stars 



N.G. Shchukina 

Main Astronomical Observatory, National Academy of Sciences, Zabolotnogo st., 27, 

03680 Kyiv, Ukraine 

O 

shchukin@mao . kiev . ua 

<N 

"o ' J. Trujillo Bueno 1 

o. 

q , Instituto de Astrofisica de Canarias, 38205 La Laguna, Tenerife, Spain 

(N 

jtb@iac . es 

' and 

o 



M. Asplund 



Research School of Astronomy and Astrophysics, Mt Stromlo Observatory, Cotter Road, 

Weston, ACT 2611, Australia 

Oh' 

i 1 martinOmso.anu.edu.au 
O 

■ ABSTRACT 



This paper presents the results of a detailed theoretical investigation of the 
impact of non-LTE effects and of granulation inhomogeneities on the derived iron 
and oxygen abundances in the metal-poor halo subgiant HD140283. Our analysis 
is based on both the 'classical' one-dimensional (ID) stellar atmosphere models 
and on the new generation of three-dimensional (3D) hydrodynamical models. 
The non-LTE calculations presented here have been carried out without inelastic 
collisions with neutral hydrogen atoms. We find that if non-LTE effects are taken 
into account when synthetizing the Fe I spectrum in both type of models, then 
the derived iron abundance turns out to be very similar in both cases. The 
emergent spectral line profiles in both models are very much weaker in non-LTE 
than in LTE because the UV overionization mechanism produces a very strong 
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underpopulation of the Fe I levels, in particular in the granular regions of the 
3D model. As a result, the non-LTE effects on the derived iron abundance are 
very important, amounting to ~0.9 dex and to ~0.6 dex in the 3D and ID cases, 
respectively. On the other hand, we find that non-LTE and 3D effects have to 
be taken into account for a reliable determination of the iron abundance from 
weak Fe II lines, because the significant overexcitation of their upper levels in 
the granular regions tend to produce emission features. As a result such Fe II 
lines are weaker than in LTE and the abundance correction amounts to ~0.4 dex 
for the 3D case. 

We derive also the oxygen-to-iron abundance ratio in the metal-poor star 
HD140283 by using the Ol triplet at 7772-5 A and the forbidden [Ol] line at 
6300 A. Our 3D results for the oxygen abundance confirm the values reported in 
some recent investigations. While the oxygen abundance derived from the O I IR 
triplet is not very sensitive to the presence of granulation inhomogeneities, such 
3D effects amount to ~ — 0.2 dex for the [Ol] line. The non-LTE abundance 
correction for the Ol IR triplet turns out to be —0.2 dex, approximately. 

Interestingly, when both non-LTE and 3D effects are taken into account there 
still remains significant discrepancies in the iron abundances derived from Fe I and 
Fe II lines, as well as in the oxygen abundances inferred from the O I and [O i] lines. 
We conclude that the discrepancies could be due to uncertainties in the stellar 
parameters of this metal poor star. We argue that adopting T c h^5600 K (instead 
of T cff ^5700 K) and [Fe/H]«-2.0 (instead of [Fe/H]«-2.5) substantially reduces 
the discrepancies in the abundances of iron and oxygen inferred from several 
spectral lines. Under such circumstances we find [O/Fe]f»0.5 at [Fe/H]=— 2. 
Obviously, our tentative conclusion that the metalicity of this type of metal- 
poor stars is significantly larger than previously thought may have far-reaching 
implications in stellar astrophysics. 

Subject headings: Galaxy: abundances — Galaxy: evolution — nucleosynthesis, 
abundances — stars: late-type stars — Sun: abundances 



1. Introduction 

Reliable determinations of the oxygen and iron abundances in the atmospheres of metal- 
poor stars are of paramount importance for 21st Century astrophysics. This is because the 
abundances of oxygen and iron in metal-poor stars represent fundamental astrophysical quan- 
tities for determining the ages of the oldest stars, for constraining the models of the chemical 
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evolution of the Galaxy, for estimating the yields of light elements (lithium, beryllium and 
boron) from spallation processes in the interstellar and circumstellar gas, for testing different 
scenarios of nucleosynthesis in supernovae, etc. Such studies largely rely upon the oxygen- 
to-iron abundance ratio ([O/Fe] = log(0/Fe) + — log(O/Fe) ) versus the stellar metallicity 
[Fe/H]. 



1.1. The oxygen abundance problem in metal-poor stars 

The shape and spread of the [O/Fe] vs. [Fe/H] relation is still hotly debated in the 
literature (e.g., Israelian et al. 2001; Asplund & Garcia Perez 2001; Melendez & Barbuy 
2002; Nissen et al. 2002; Fulbright & Johnson 2003; Takeda 2003, and references therein). 

The controversy is linked with discordant results obtained using different indicators of 
the oxygen abundance, such as: i) the oxygen [O i] forbidden A6300 A line; ii) the oxygen O I 
infrared triplet at 7772-5 A; iii) the OH lines in the near ultraviolet; iv) the OH lines in the 
infrared. More precisely, the oxygen abundance derived using all these lines agree more or less 
with each other for moderate metal-deficient disk stars having —1 < [Fe/H] < 0. However, 
there are clear discrepancies among the results obtained for metal-poor halo stars with lower 
metallicities. Some authors find that the [O/Fe] ratio shows a plateau between 0.4 and 0.5 
dex in the metallicity range -3<[Fe/H]< - 1 (Barbuy 1988; Kraft et al. 1992; Fulbright & 
Kraft 1999; Asplund & Garcia Perez 2001; Caretta et al. 2000; Melendez & Barbuy 2002), 
while others find a linear increase with decreasing metallicity, reaching [0/Fe]~ +1 dex 
at [Fe/H]= -3 (Abia & Rebolo 1989; Tomkin et al. 1992; Cavallo et al. 1997; Israelian 
et al. 1998, 2001; Boesgaard et al. 1999; Mishenina et al. 2000; and references therein). 
Recently, Nissen et al. (2002) have analyzed [O i] and [Fe n] lines within the framework of 
3D atmospheric models and LTE line formation concluding that there is a quasi-linear trend 
for [Fe/H]< -2, with [0/Fe]~ +0.5 dex at [Fe/H]= -2.5. Furthermore, an extensive ID 
re-analysis of published equivalent widths for the O I IR and the forbidden [O i] lines (see 
Takeda 2003) shows that the discrepancy between O I and [O i] lines in metal-poor halo stars 
tends to be smaller for stars with higher effective temperatures (T cff ) and gravities (logg), i.e. 
preferentially for dwarfs and subgiants than for giants. Fulbright & Johnson (2003) support 
this conclusion. In addition, the oxygen abundances in dwarfs derived from the O I IR lines 
tend to have values between those obtained with the UV OH and [O i] lines. 
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1.2. The iron abundance problem in metal-poor stars 

A key problem directly related with the [O/Fe] issue in metal-poor stars is the reliability 
of the iron abundance determinations. In fact, the derivation of the iron abundance from Fe I 
and/or Fell lines is a very complex problem that requires detailed non-LTE investigations 
(hereafter NLTE) like the one carried by Shchukina & Trujillo Bueno (2001) and in this 
paper. In general, the Fe I lines suffer from UV overionization in the atmospheres of late- 
type stars (Athay and Lites 1972; Rutten 1988; Shchukina & Trujillo Bueno 2001). This 
NLTE effect tends to be substantially more important in metal-poor stars than in solar-like 
stars as a result of the lower electron density and of the weaker UV blanketing (Thevenin 
& Idiart 1999; Shchukina and Trujillo Bueno 2001; Israelian et al. 2001, 2004; Shchukina et 
al. 2003). We should mention also that Gratton et al.(1999) concluded that NLTE effects 
are rather small in metal-poor stars. However, their calculations presumably overestimated 
the efficiency of inelastic collisions with hydrogen atoms which thermalize the Fe I levels and 
thus almost restore LTE conditions. 

It is of interest to point out that while the Fel lines are (in principle) expected to 
be sensitive to the UV overionization mechanism, the LTE approximation is thought to be 
suitable for weak lines of Fell (see, e.g., Asplund et al 1999; Nissen et al. 2002). However, 
as we shall see below this does not have to be necessarily the case because many subordinate 
Fe II lines whose upper levels are of odd parity (starting at ^ 6 D°) might be partially filled by 
emission due to an optical pumping mechanism similar to that investigated by Cram et al. 
(1980) for the solar case. 

1.3. The solar oxygen abundance problem 

Obviously, the solar case constitutes a fundamental reference for a reliable quantification 
of chemical abundances in other stars. The solar iron abundance appears to be fairly well 
determined now (loge Q (Fe)=7.50±0.10; see Shchukina & Trujillo Bueno 2001) — that is, it 
is similar to the meteoritic value. In fact, this is the value adopted by Nissen et al. (2002) 
in their study of the O/Fe ratio in metal-poor stars. 

Great progress has been also made towards a reliable determination of the solar oxygen 
abundance. The ID modeling approach leads to discrepancies between the solar oxygen abun- 
dance determined from the permitted O I lines and from the [O i] and OH lines. Abundance 
determinations from OH lines using the ID model of Holweger & Miiller (1974) suggested 
first a high oxygen abundance value (loge (O)=8.93±O.O3; see Grevesse et al. 1984; Sauval 
et al. 1984; Anders & Grevesse 1989). Several years later, Grevesse & Sauval (1998) obtained 
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an oxygen abundance that was lower by 0.1 dex. This lower value was found after using im- 
proved loggf values and a slightly modified ID model, which helped to remove the existence 
of a dependency of the derived oxygen abundance with the excitation potential of the lower 
level of the OH transitions. The forbidden [O i] line gives also a high oxygen abundance 
value when the contribution of the blend due to a Nil line is ignored (Lambert, 1978; Sauval 
et al. 1984; Holweger 2001). However, the permitted Ol lines at A7772-5 A yield a value 
for the oxygen abundance that is lower by 0.2 dex (Kiselman 1993; Holweger 2001). Some 
authors have argued that the existing NLTE calculations for these spectroscopic indicators 
of the oxygen abundance might not be valid due to the neglect of inelastic collisions with 
hydrogen. For this reason the solar oxygen abundance derived from the [O i] forbidden line 
was considered as the most reliable one, because of the suitability of the LTE approximation 
(hence, uncertainties of collisional rates with hydrogen atoms do not play any role). Another 
argument used to support a high value for the oxygen abundance is the small scatter in the 
oxygen abundances obtained from different OH lines. Kiselman & Nordlund (1995) tried 
to clarify the solar oxygen abundance issue by using a 3D solar photospheric model that 
is similar but less sophisticated than that of Asplund et al. (2000). They concluded that 
a downward revision of the solar oxygen abundance is required, with respect to the value 
derived from [Ol] and OH lines. More recently, Allende Prieto et al. (2001) have obtained a 
more precise estimate by taking into account the influence of the above-mentioned blend by 
the Nil line and by using the 3D hydrodynamical model of Asplund's et al. (2000). They 
found loge (O)=8.69±O.O5, which shows that the [Ol]-based abundance derived via the 3D 
approach is significantly lower than that derived from classical ID model atmospheres. 

In a recent IAU Symposium on modeling of stellar atmospheres Shchukina et al. (2003) 
reported that the oxygen abundance obtained via NLTE synthesis of the solar O I IR triplet 
in the 3D solar photosphere model of Asplund's et al. (2000) is loge (O)~8.7. This result is 
in good agreement with both Allende Prieto's et al. (2001) conclusion and with the recent 
3D analysis of Asplund et al. (2004) who found excellent agreement between the averaged 
oxygen abundance values derived from [O i] , O I, OH vibration-rotation and OH pure rotation 
lines 1 . 



1 It is interesting to mention that the oxygen abundances obtained by Asplund et al. (2004) from OH 
rotational lines using the 3D modeling approach show a large line-to-line scatter, and with sensitivity to the 
line strength and to the excitation potential of the transition's lower level. However, such presumably minor 
discrepancies are not the topic of this paper. Here we restrict ourselves to consider the NLTE determination 
of the solar oxygen abundance via O I lines, for consistency with similar computations we have carried out 
for the metal-poor star HD140283. 
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1.4. Possible reasons for the abundance discrepancies 

It is currently believed that abundance discrepancies in metal-poor halo stars and the 
sun are dominated by deficiencies in the spectral line modeling. Errors in the measurements 
of equivalent widths are of minor importance with the exception of weak [O i] and Fe II 
lines. All abundance indicators in metal-poor stars are to a different degree sensitive to the 
adopted stellar parameters (T eS , \ogg and [Fe/H]). They have to depend also on atmospheric 
properties such as the temperature gradient and the ratio of the line to continuous opacity. 
Large uncertainties could result from neglecting atmospheric inhomogeneities and velocity 
fields (i.e., granulation inhomogeneities produced by the stellar surface convection process). 
Part of the discrepancies are probably due to NLTE effects. There is now an abundance of 
literature concerning the iron and oxygen line formation problems in classical plane-parallel 
stellar atmospheric models (Thevenin & Idiart 1999; Mishenina et al. 2000; Kiselman 2001; 
Takeda 2003; and references therein). The recent generation of 3D hydrodynamical models 
of stellar atmospheres have activated the present debates on the influence of granulation 
inhomogeneities and of NLTE effects on the determination of chemical abundances in metal- 
poor stars (see Asplund et al. 1999; Asplund and Garcia Perez et al. 2001; Shchukina & 
Trujillo Bueno 2001; Nissen et al. 2002; Asplund et al. 2004). In this respect, it is very 
important to emphasize that the 3D hydrodynamical simulations of Asplund and coworkers 
show that the atmospheres of metal-poor halo stars may have substantially lower kinetic 
temperatures than the radiative equilibrium temperatures of the 'classical' theoretical ID 
model atmospheres of Kurucz (1993) and the MARCS collaboration (Gustafsson et al. 1975). 



1.5. Motivation 

The previous discussion makes evident that detailed, 3D+NLTE investigations of the 
iron and oxygen abundances are urgently needed. However, it is crucial that this type of 
analyses are carried out jointly for iron and oxygen, and using the same set of atmospheric 
models, continuum opacity package and radiative transfer code. The present paper presents 
the results of such an investigation. We study the NLTE formation problem of the Fel, 
Fe II and O I lines in a 3D hydrodynamical model of the atmosphere of the metal-poor star 
HD 140283. This classical halo subgiant is widely used for studies of the very early Galaxy 
and the origin of the chemical elements. In order to quantify the impact of 3D simulations 
on the abundance determinations in HD140283 we carry out a detailed comparison with 
the ID classical modeling. To this end, we use the 3D hydrodynamical atmospheric model 
of Asplund's et al. (1999) for the star HD140283, as well as a suitable grid of Kurucz's 
ID models. Moreover, we evaluate the extent to which 3D+NLTE effects influence the 
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determination of the stellar parameters of this metal-poor subgiant. 

Our paper considers also the case of the Sun as a reference for stellar abundance deter- 
minations with emphasis on the solar oxygen abundance derived from the O I IR triplet. As 
in our NLTE+3D analysis of the solar iron abundance (Shchukina & Trujillo Bueno 2001) we 
investigate here the solar oxygen abundance problem by using the same 3D hydrodynamical 
solar granulation model. 

2. The diagnostic method 

In order to determine the stellar abundances of iron and oxygen we have to compare 
some of the properties of the observed spectral line profiles (e.g., their equivalent widths) 
with those calculated via the numerical solution of the radiative transfer equation in a given 
model of the atmosphere of the star under investigation. To this end, we have used the 
same NLTE code and strategy that we applied in our previous investigation of the iron line 
formation problem in a 3D hydrodynamic model of the solar photosphere (see Shchukina & 
Trujillo Bueno 2001), with the only difference that we have now computed flux profiles. Our 
3D strategy for investigating the issue of the abundances of iron and oxygen in the metal- 
poor star HD140283 is similar to that of described by Shchukina & Trujillo Bueno (2001) 
for solar-like atmospheres, but using a 3D snapshot model from Asplund's et al. (1999) 
stellar surface convection simulation. Therefore, we have neglected the effects of horizontal 
radiative transfer on the atomic level populations -that is, we have used the so-called 1.5D 
approximation. As discussed by Shchukina & Trujillo Bueno (2001) this is expected to be 
a good approximation for determining the stellar iron abundance from Fel and Fell lines. 
Kiselman & Nordlund (1995) already demonstrated that it is also a good approximation for 
the oxygen IR triplet. 

Our NLTE code is based on very efficient iterative methods (Trujillo Bueno & Fabiani 
Bendicho 1995; Socas-Navarro & Trujillo Bueno 1997; Trujillo Bueno 2003) that were de- 
veloped to facilitate radiative transfer simulations with complex atomic systems in realistic 
stellar atmospheric models. We have carried out both LTE and NLTE spectral line syntheses 
by using both the ID and 3D models mentioned below. In each case, the abundances of iron 
and oxygen are modified iteratively till achieving the best fit to the observations. 
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3. The observations 

We have chosen 33 Fe I and 15 Fe II lines for our 3D and ID determinations of the iron 
abundance in the halo subgiant HD140283. The equivalent widths {W) of all the Fe I lines and 
of 11 Fell lines were obtained from high- resolution spectra observed a few years ago with 
the Harlan J. Smith 2.7m telescope at McDonald observatory (R. Garcia Lopez, private 
communication; see also Allende Prieto et al. 1999). The equivalent widths (W) for the 
remaining 4 Fe II lines are taken from Table 3 of Nissen et al. (2002). They are weak lines with 
W in the range 1-4 mA . Their wavelengths are AA6149.238, 6247.562, 6432.654, 6456.376 
A. The iron lines used in our analysis are listed in Table 1. 

The Ol IR triplet lines in the metal-poor star HD 140283 are weak. Given uncertainties 
in the measurements we have averaged their equivalent widths over three sets of observations 
(Abia & Rebolo 1989; Tomkin et al. 1992, and G. Israelian, private communication). The 
resulting equivalent widths of the 7772, 7774 and 7775 A lines used in our oxygen abundance 
determinations are 7.9, 4.8, 3.4 mA, respectively. The equivalent width of the forbidden [O i] 
A6300 A line (W = 0.5 mA) has been taken from Nissen's et al. (2002) paper. The oxygen 
line data are specified in Table 1. 

Measurements of the equivalent widths of the solar O I IR lines tend to have an un- 
certainly of ~5%. Therefore, we have used averaged values of five sets of observations at 
the solar disk centre (Jungfraujoch atlas of Delbouille et al. 1973; Gurtovenko & Kostik 
1989; Altrock 1968; Miiller et al. 1968; King & Boesgaard 1995). The resulting averaged 
solar equivalent widths for the 7771.96, 7774.18 and 7775.40 A lines are 81.6, 70.4, 57.2 mA, 
respectively. 

4. The atmospheric models 

In each case (i.e., the solar and the metal-poor star) we have used a single 3D snapshot 
model taken from realistic radiation hydrodynamical simulations of stellar surface convection 
(Stein & Nordlund 1998; Asplund et al. 1999, 2000). Such time-dependent simulations 
have no adjustable free parameters, since they are based on the numerical solution of the 
equations of mass, momentum and energy conservation for given values of the effective 
temperature (T e g), surface gravitational acceleration (logg) and chemical composition. In 
order to facilitate the radiative transfer calculations the original snapshots, which have 82 
gridpoints along the vertical direction, were interpolated to grids of 50 x 50 x 102 and 
50 x 50 x 121 points for the solar and metal-poor star cases, respectively. The improvement 
in the vertical resolution (i.e., with 102 and 121 gridpoints sampling the line and continuum 
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forming regions) was introduced to enhance the accuracy of the NLTE calculations. 

The solar 3D model is that described by Asplund et al. (2000). This 3D model of 
the solar photosphere has been used by Shchukina & Trujillo Bueno (2001) in their NLTE 
investigation of the solar iron abundance. We have also used it to determine the solar oxygen 
abundance via NLTE analysis of the O I IR triplet at 7772-5 A. 

Information on the corresponding 3D model for the star HD140283 can be found in 
Asplund et al. (1999). The adopted stellar parameters were logg = 3.7, T cS = 5700 K and 
[Fe/H]=-2.5, which result from the infrared flux method (IRFM) for T eS , from Hipparcos 
parallaxes for logg and from published values for the [Fe/H] (see Asplund et al. 1999). In 
our differential comparison between the results obtained via the ID and 3D approaches we 
have used atmospheric models characterized by the same stellar parameters (T e g, logg and 
[Fe/H]). In Sect. 8 we will also show results in a grid of 'classical' ID models spanning the 
range 5300 < T cS < 5700, 3.7 < hgg < 3.9, -2.0 < [Fe/H] = -2.5. 

5. The atomic models 

In our spectral synthesis of oxygen and iron lines we adopted the quantum-mechanical 
approach developed by Anstee & O'Mara (1995); Barklem & O'Mara (1997), and Barklem 
et al (1998) for the Van der Waals collisional broadening. In the following subsections we 
provide some information about the atomic models of oxygen and iron we have used for our 
radiative transfer calculations. 



5.1. Oxygen 

Our atomic model for oxygen is based on the data of Carlsson & Judge (1993). It has 23 
O I fine structure levels and one O II level. We considered 31 bound-bound and 23 bound- free 
radiative transitions. The O I Grotrian diagram has a helium-like structure. The triplets and 
the quintets are only weakly coupled to each other. The two systems differ in one important 
aspect: the triplets are connected to the ground state by strong resonance lines, while the 
quintets are not. The latter form a separate system with the 3 5 S level acting as a ground 
level and the 7772-5 A IR triplet as resonance lines. O I levels have hydrogen-like energies. 
As a result, there is an important wavelength coincidence of the Ol 1025 A line with the 
hydrogen Lyman — j3 line. We emphasize that the NLTE results for O I are quite insensitive 
to including additional levels and transitions (Shchukina 1987; Kiselman 1993; Takeda 2003). 
In agreement with previously published results we have found that the effects of the coupling 
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between the triplet and quintet systems, Ly^ pumping, and the binding of CO molecules 
are marginal in the stellar atmospheres studied here (see Shchukina 1987; Kiselman 1993; 
Carlsson & Judge 1993; Kiselman & Nordlund 1995). 

5.2. Iron 

Our model atom for Fe I+Fe Il+Fe ill is similar to that used by Shchukina & Trujillo 
Bueno (2001). It has 225 Fel levels, 23 Fell levels and 1 level for Fein including their 
multiplet fine structure. The Fel levels are interconnected by 330 bound-bound radiative 
transitions, while we have 25 strong radiative transitions among the Fe II levels. All levels are 
coupled via collisional transitions with electrons. The Fel term diagram is, in fact, complete 
up to an excitation potential x = 5-72 eV. At higher energies it contains about 50% of the 
terms that have been presently identified. Each term of our atomic model is coupled to its 
parent term of the next ionization stage by photoionization transitions. For most of the Fe I 
terms that have x <5 eV we used the experimental photoionization cross sections given by 
Athay & Lites (1972). For the rest of them, and also for the Fell terms, the hydrogen-like 
approximation was used. 2 

In our calculations we have not included inelastic collisions with hydrogen atoms. We 
have also neglected the UV haze opacities. We have carried out several numerical exper- 
iments in order to investigate how large could be the uncertainties in the iron abundance 
determinations due to uncertainties in the hydrogen collisional rates, in the UV haze opacity 
and in the bound-free cross sections. Considering all factors together we have found that 
the uncertaintly for the solar iron abundance does not exceed 0.07 dex (see Shchukina & 
Trujillo Bueno 2001). The problem of the iron abundance uncertainty for the metal-poor 
star HD 140283 is discussed in Section 8. 



6. The oxygen abundances 

In this section we report on the results we have obtained for the oxygen abundance in 
the sun and in the metal-poor star HD140283. Our abundance determinations and some 
recent abundance studies by other authors are briefly summarized in Tables 2 and 3. 



2 It is interesting to mention that the new radiative cross-sections of Bautista (1997) are substantially 
larger than those given by the hydrogenic approximation. Therefore, in what respects the influence of 
photoionization cross sections, we might perhaps be underestimating the NLTE effects. 



6.1. The Sun 



We have derived the abundance of oxygen in the solar photosphere by fitting the ob- 
served disk centre profiles and equivalent widths of the O I IR triplet via spectral synthe- 
sis in the 3D model of the solar photosphere. Taking into account NLTE effects we find 
loge(O) = 8.70 ± 0.06 (see also Shchukina et al. 2003), while the LTE approximation gives 
loge(O) = 8.93 ± 0.06. Note, that the uncertainties reflect the line-to-line scatter for the OI 
IR triplet lines. In both cases we obtained an almost perfect fit to the observed disk centre 
profiles and equivalent widths. Our NLTE result is in good agreement with that obtained by 
Allende Prieto et al. (2001) from the forbidden [Ol] 6300 A line (loge(O) = 8.69). It is also 
consistent with the 3D result of Asplund at al. (2004) for the OH lines (loge(O) = 8.61 — 8.65) 
and with their full 3D+NLTE result for the Ol triplet (loge(O) = 8.65). It is important to 
note that the small discrepancy between ours and Asplund's et al. (2004) solar abundance 
value derived from the O I IR triplet is mainly due to differences in the oscillator strengths 
used (Alogg/ —0.05 dex). The oxygen abundance changes that result from a 5% un- 
certainty in the solar equivalent width measurements of these spectral lines lie in the same 
range (ss ±0.05 dex). 

The NLTE formation of the O I triplet can be easily understood via a two-level atom 
approach if the line opacity is calculated assuming LTE (Shchukina 1987; Kiselman 1993; 
Kiselman & Nordlund 1995; Takeda 2003; Shchukina et al. 2003, Asplund et al. 2004, 
and references therein). The photon losses mechanism lowers the line source functions both 
in the granular and intergranular regions of the 3D model, which implies that the NLTE 
profiles of the IR triplet lines are deeper than in LTE. The effect is larger in the inter- 
granular regions, but the spatially averaged emergent profiles are only slightly affected by 
the granulation inhomogeneities. Consequently, the effects of granulation inhomogeneities 
on the NLTE solar oxygen abundance determination are small, i.e. ~ 0.06 dex (Shchukina 
et al. 2003). More significant are the NLTE abundance corrections for the Ol infrared 
triplet lines observed at the solar disk centre, which amount to —0.25, —0.22, —0.20 dex 
for AA7771. 96, 7774.18, 7775.40 A, respectively — that is, the LTE approximation tends to 
overestimate the derived oxygen abundance. Such NLTE effects are very similar to those 
reported by Asplund et al. (2004), which result from full 3D+NLTE computations. 



6.2. The metal poor halo subgiant HD140283 

It is convenient to begin this section by giving the values of the oxygen abundances 
we have obtained without assuming LTE when using the ID and 3D models of this metal 
poor star. As mentioned above, both models are characterized by \ogg = 3.7, T cff = 5700 K 
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and [Fe/H]= -2.5. The ID model gives loge(O) = 7.11 ± 0.06 for the Ol IR triplet, while 
loge(O) = 7.18 for the [O i] forbidden line at 6300 A. Concerning the 3D-based determination 
we find loge(O) = 7.08 ± 0.06 when using the Ol IR triplet, but loge(O) = 6.96 for the 
[Ol] line. Our ID results differ slightly from those presented by Nissen et al. (2002), 
probably because we have not used the same continuum opacity package and/or identical 
ID atmospheric models. Our estimates indicate that the abundance differences could reach 
~ 0.05 - 0.07 dex. 

Having given these values, it is now of interest to note that the NLTE effects in the ID 
and 3D atmospheric models of this metal poor star are not negligible, but smaller than in 
the solar case because the O I IR lines originate in deeper layers where their source functions 
and the level populations of their upper levels are closer to their LTE values (Shchukina et 
al. 2003). Adopting LTE for the Ol lines leads to an abundance overestimation of 0.18 and 
0.16 dex for the ID and 3D cases, respectively. The determination of the oxygen abundance 
based on the O I IR triplet is nearly insensitive to 3D effects. The correction does not exceed 
0.03 dex. 

On the contrary, the forbidden [O i] line at 6300 A is sensitive to the stellar atmospheric 
inhomogeneities. In agreement with Nissen et al. (2002) we find that a 3D analysis reduces 
the oxygen abundance derived from this forbidden line by 0.22 dex. The NLTE correction 
is extremely small (< —0.01 dex). 

7. The iron abundance 

This section focusses on the iron abundance in the Sun and in HD 140283. Our results 
and a comparison with those of recent publications are given in Tables 2 and 3. 

7.1. The Sun 

Shchukina & Trujillo Bueno (2001) applied a 3D+NLTE approach to determine the 
solar iron abundance from 18 Fel lines that were carefully selected. Ideally, one would like 
to determine the iron abundance in the Sun and HD140283 employing the same Fel line list. 
However, due to the relatively high solar metallicity the solar spectral lines at wavelengths 
below 5000 A are much more affected by blends than in the metal-poor star HD 140283. As 
a result, overlapping wings of nearby strong lines cause large uncertainties in the equivalent 
width measurements of most of the solar Fel lines listed in Table 1. The only exception are 
the lines at A > 5700 A (i.e. only 8 Fel lines from 33 ones). In any case, the main reason 
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why not the same lines can be used in the Sun and HD140283 is that the line strengths are 
vastly different in the two cases. Therefore, the weak iron lines in HD 140283 are very strong 
and therefore ill-suited for abundance purposes in the Sun. Thus, to recalculate the iron 
abundance of the Sun using the lines of Table 1 is hardly possible. 

The solar Fel lines used in our 3D solar abundance determination (see Table 1 in 
Shchukina & Trujillo Bueno 2001) were selected from the well-known "Kiel+ Oxford" line 
list according to the criteria discussed by Kostik, Shchukina & Rutten 1996. The lines chosen 
by Shchukina & Trujillo Bueno (2001) account for the iron abundance behaviour versus the 
lower-level excitation potential x -that is, for the existence of both an abundance dependence 
on x an d a sizable scatter in the abundance values for each x- Such 18 representative Fe I 
lines yield the largest and lowest abundances for each %-value and spaning a large range in 
X which goes from to 4.6 eV. This should be a very suitable choice of lines in order to 
evaluate the largest 3D effects. Shchukina & Trujillo Bueno (2001) obtained a Fe abundance 
loge Q (Fe)=7.50±0.10 for the 3D+NLTE case. The 3D+LTE approach gives a lower abun- 
dance, i.e. loge Q (Fe)=7.43±0.11. Neglecting the effect of granulation inhomogeneities leads 
to a considerably larger value (loge©(Fe)=7.62±0.04; e.g. Kostik et al. 1996). 

We stress that the Fel lines we have used for the star HD140283 were chosen to obtain 
a combination of lines that is comparable to that provided by the 18 solar lines (concerning 
their equivalent widths and x coverage). 

7.2. The metal-poor halo subgiant HD140283 

Like in the solar case (Shchukina and Trujillo Bueno 2001), near-UV overionization 
appears to be the main NLTE mechanism for Fel in this metal-poor star. Both in the 
ID and 3D models it causes an important underpopulation of the Fel levels and hence a 
reduction in the opacity of the Fel lines (Fig. 1, left panel). As a result, the atmospheric 
heights where the line optical depth reaches unity are shifted toward deeper layers where 
the line source functions turn out to be close to the Planck function. The relevance of this 
effect increases with the vertical temperature gradient, which is larger above the granules. 
For this reason, we find a strong weakening of the Fe I lines that emerge from such granular 
regions, while above the intergranular lanes the emergent spectral line profiles do not change 
appreciably. On the whole, the calculated NLTE profiles turn out to be much weaker than 
when the LTE approximation is used (see Fig. 2). 

Concerning NLTE effects in Fell lines it is very important to point out that in the 
granular regions we find a significant overpopulation of the excited levels of odd parity 
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starting at z 6 B° (Fig. 1, right panel). According to Cram et al. (1980) and Rutten (1988) 
this effect is caused by an excess in the UV continuum radiation field near 2600 A, which leads 
to a pumping of the upper level populations of the Fe II lines. As a consequence, the source 
functions of the Fe II lines with high excitation energy of the upper level (x U p > 4.8 ev) exceed 
significantly the Planck function values. This implies that such Fell lines are weaker than 
stipulated by the LTE approximation and may even be in emission instead of in absorption. 
This radiative pumping mechanism takes place in both the ID and 3D models. Like in the 
case of Fe I lines it is sensitive to the vertical temperature gradient. The steeper the gradient, 
the larger the pumping. As a result, in the granular regions where the temperature gradients 
are larger the effect is much more pronounced (see Fig. 1). 

The importance of such an effect is illustrated in Fig. 3, which shows both NLTE and 
LTE results for two Fell lines using the 3D model atmosphere of the metal-poor subgiant 
HD 140283. We point out that emission is much more pronounced for the emergent pro- 
files that originate above the stellar granular regions, while above the intergranular lanes 
the spectral line profiles are hardly modified with respect to the LTE case. Allowing for 
departures from LTE produces a drastic change in the calculated spatially averaged flux of 
the weak Fell A6247.56 A line in the 3D model: the absorption profile is replaced by one in 
emission. Note, that this spectral line is taken from the set of weak Fell lines that Nissen 
et al. (2002) assumed to be insensitive to NLTE effects. However, we find that such type 
of lines do suffer from the above-mentioned UV pumping mechanism because all of them 
belong to the aforementioned category of lines, i.e. their upper levels are of odd parity and 
are situated above the z 6 D°-level. The stronger Fell A5316.62 A line displayed in Fig. 3 is 
also from this category. Figure 3 clearly demonstrates that the stronger the spectral line the 
smaller the relative effect of the UV-pumping. For this line emission arises only at the core 
of the "granular" profiles. Spatial averaging masks completely such an emission, but makes 
the flux profile weaker than in LTE. It is also of interest to mention that in the ID stellar 
atmosphere model departures from LTE produce similar changes in the Fel and Fell lines 
but they are less pronounced than in the 3D case. For instance, the flux profile of the weak 
Fell A6247.56 A line turns out to be in absorption for the 1D+NLTE case. 

How large could be the impact of UV-overionization and UV-overexcitation on the iron 
abundance determination in the star considered? It is obvious that in order to compensate 
for a weakening of the emergent spectral line profiles in the stellar atmosphere model un- 
der consideration the iron abundance must be increased in comparison with the LTE case. 
Figure 4 quantifies this conclusion. It shows the LTE versus NLTE iron abundances derived 
from either Fe I or Fe II lines in the 3D and ID stellar atmospheric models of the metal poor 
star HD140283. As seen in the figure, the NLTE abundance corrections for the Fel lines 
are extremely large in the 3D model (~ 0.9 dex) and a bit smaller in the ID case (~ 0.6 
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dex). On the other hand, it is very important to point out that NLTE effects on the iron 
abundance are most pronounced for the weak Fe II lines. The four lines used by Nissen et al 
(2002) for their iron abundance determinations reveal strikingly large deviations from LTE. 
On average, the NLTE abundance corrections for the Fell lines increase from ~ 0.16 dex in 
the ID case to ~ 0.4 dex in the 3D model. 

When NLTE effects are taken into account we find no difference between the mean 
Fel-based abundances obtained via the ID and 3D models (loge(Fe)~5.77). Therefore, 
when NLTE effects are taken into account the Fel lines formed in the 3D model are hardly 
sensitive to granulation inhomogeneities. However, for the Fell lines the 3D abundance 
correction (3D+NLTE - 1D+NLTE) is significant — that is, 0.26 dex. 

When the LTE approximation is used the mean iron abundance obtained via spectrum 
synthesis of Fel lines in the 3D model is lower than for the ID case (cf. Asplund et al. 1999; 
Nissen et al. 2002). The differences between the ID and 3D determinations depend on the 
lower excitation potential x of the spectral line under consideration. They range from ~ 0.6 
dex at x — eV to ~ 0.2 dex at x ~ 2 eV, while the mean difference is smaller (~ 0.3 dex). 
In contrast, on average the Fell lines are insensitive to 3D effects when LTE is assumed. 

8. Revising the stellar parameters 

There are several discrepancies indicating that the iron abundance obtained in Sect. 7.2 
for the star HD 140283 cannot be considered as reliable, neither in the ID case nor in the 3D 
case. For example, we have (1) a divergence of the average abundances derived from both 
ionization stages, (2) a correlation of the abundances derived from Fe I lines with the lower- 
level excitation potential and (3) a large scatter in the abundances derived from Fell lines 
(particularly, the weak ones). The elimination of any of such discrepancies cannot be simply 
achieved in terms of uncertainties in observed equivalent widths, oscillator strengths, inelastic 
electron collisional rates and photoionization cross sections. Our numerical experiments 
show that the typical errors introduced by uncertainties in such quantities do not exceed 0.1 
dex (see Shchukina and Trujillo Bueno 2001). The iron abundance error that results from 
neglecting the UV haze opacity for the metal-poor star HD140283 is even smaller than 0.1 
dex. In our opinion, a reasonable possibility for explaining the reported discrepancies is to 
revise the stellar parameters of the subgiant HD140283. 

In this section we aim at finding the "best-choice" of stellar parameters for the star 
HD 140283, both for the NLTE and LTE cases. It is important to emphasize from the 
outset that in our NLTE computations we have neglected inelastic collisions with neutral 
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hydrogen atoms for excitation and ionization. At present, there are no reliable values for 
the ensuing collisional rates. The often-used classical Drawin (1968) formula leads to very 
uncertain estimates and to a substantial overestimation of the collisional rates when detailed 
quantum mechanical and/or laboratory data is available, such as for Li and Na (see Biberman 
et al. 1987; Lambert 1993; Holweger 1996; Barklem et al. 1998, 2003; Belyayev et al. 
1999; Shchukina & Trujillo Bueno 2001; Asplund et al. 2004), which explains why the 
conclusions of Gratton et al. (1999) have been critized. Obviously, inelastic collisions with 
neutral hydrogen atoms would tend to reduce the NLTE effects. Therefore, the difference 
between the NLTE and LTE iron abundances may be considered as the maximum effect 
that inelastic collisions might produce. In our opinion, the tentative conclusion of Korn, 
Shi & Gehren (2003), achieved via a 1D+NLTE analysis, that hydrogen collisions might 
be efficient in the atmosphere of HD 140283 needs to be re-investigated in the light of 3D 
hydrodynamical models because, as pointed out by these authors, the lower temperatures of 
the 3D models would alleviate their need of invoking thermalizing collisions to counterbalance 
photoionization. 

Our analysis for the revision of the stellar parameter is based on the following criteria: 

i) the average abundances obtained from Fe I and Fe II lines have to be equal. 

ii) the "best choice" solution has to give the minimum mean standard deviation (a) for 
the abundances derived from Fel and Fell lines. 

iii) abundances have to be independent of the lower excitation potential (x) ■ We restrict 
ourselves to revise only the metallicity and the effective temperature. The given surface 
gravity of the star HD 140283 is reliable, since it has been obtained from accurate Hipparcos 
parallaxes. 

Figure 5 helps us to find the "best-choice" of stellar parameters for the star HD 140283. 
The left panels show the variation of the Fel and Fell abundances with metalicity, T ef f, 
and logg in a suitable grid of ID model atmospheres. Our results clearly demonstrate that 
changes in T eff seriously affect the Fe i-based abundances. However, such changes in T cS are 
unimportant for abundance determinations from Fell lines. On the other hand, the Fell- 
based abundances are slightly sensitive to the changes in logg, which is not the case for the 
abundances obtained from Fell lines. In addition, the dependence on metallicity is more 
pronounced for Fel lines, while for Fell lines there is no such a variation. It is interesting 
to note that when LTE is assumed, then only the correlation of the inferred Fe I abundances 
with T eff is present. 
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8.1. Stellar parameter determinations: ID case 

Consider first the 1D+LTE modeling case for the iron lines. Figure 5 (middle panel to 
the Ihs) and Fig. 4 (bottom panel to the rhs) show that in this case identical abundances 
from both Fel and Fell lines could be obtained around T c g close to 5700 K. In fact, our 
detailed analysis gives T eff ~ 5670 K and [Fe/H]~ —2.3. With these stellar parameters Fell- 
based abundances turn out to be insignificantly scattered (a ~ 0.1 dex) around the average 
value. However, the abundances derived from Fe I lines show a sensitivity to the lower-level 
excitation potential (x) that is similar to the trend seen in the bottom panel to the rhs of 
Fig. 4. 

Secondly, consider the 1D+NLTE modeling case. In order to reach a good agreement 
between NLTE abundances derived from Fel and Fell lines we would have to reduce the 
effective temperature down to T e s = 5450 K and to increase the stellar metallicity up to 
[Fe/H]~ —2.1. With such stellar parameters the x-dependence of the Fel abundances van- 
ishes. However, the abundances derived from Fell lines become undesirably spread around 
the mean value (a ~ 0.3 dex). 

8.2. Stellar parameter determinations: 3D case 

In order to achieve a reasonable estimation of the "best choice" of stellar parameters 
the curves of Fig. 5 (left panels) have to be corrected for the 3D effects of granulation 
inhomogeneities. This means that we need to know how the 3D abundance corrections change 
with [Fe/H], T c g and logg for both the LTE and NLTE cases. Unfortunately, a suitable grid 
of 3D models for several effective temperatures and metalicities similar to that provided by 
Kurucz for the ID case is not yet available, besides the fact that NLTE computations in 
3D are presently very computationally costly. However, we can try to estimate the true 
stellar parameters of the star HD140283 via the following strategy, which is inspired by our 
conclusion that 3D effects are important for the Fe I lines only if LTE is assumed and that 
3D effects are relevant for the Fe II lines only when we allow for departures from LTE. 

Firstly, let us choose the 3D+LTE modeling approach. As we have seen, if LTE is 
assumed then the Fe I lines are sensitive to the 3D effects while the Fe II lines are not. 
Therefore, if the LTE approximation is used it is much safer to revise the stellar parameters 
via iron abundances derived from Fe II lines. We find T e g R3 5700 and [Fe/H]» —2.3 (see the 
bottom panel to the rhs of Fig. 4). 

Secondly, consider the 3D+NLTE approach. It is obvious that in this case we better 
avoid the use of Fe II lines because we have shown that they are sensitive to both 3D 
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and NLTE effects. However, we have demonstrated also that at least for the available 
3D model of the star HD 140283 the NLTE iron abundance derived from Fe I lines is very 
similar to that obtained via the 1D+NLTE modeling approach. For this reason, we think 
that we can estimate the true stellar parameters by choosing Fe I lines and a 1D+NLTE 
modeling approach. We find that the standard deviation o minimizes at T c g ~ 5600 K 
and [Fe/H]~ —2. Now the mean iron abundance turns out to be loge(Fe) = 5.54 ± 0.14. 
Moreover, the dependence of the derived iron abundances with x becomes negligible. 

8.3. Uncertainties in the determination of the stellar parameters and the 

oxygen abundance 

As we have just seen, depending on the approach used for modeling the iron lines we 
get four "best choice" sets for the stellar parameters of the star HD140283, with T e g- in the 
range ~ 5450 - 5700 K and [Fe/H] between -2.0 and -2.3 (see Table 4). How large could 
be the impact of such uncertainties in the stellar parameters of the star HD140283 on the 
determination of the oxygen abundance? The right panels of Fig. 5 clarify this issue. We 
can see that when ID models are used then both the Ol and [Ol]-based abundances are 
sensitive to changes in T e g, but not too much to changes in logg and/or in metallicity. The 
sensitivity of the NLTE abundance corrections to these parameters (logg and metallicity) 
is much smaller. On average, they vary between 0.17 and 0.19 dex for Ol lines, while it is 
negligible for the [Ol] line (< 0.01 dex). Our 1D+NLTE modeling of oxygen lines shows 
that abundances derived from Ol and [Ol] lines give the same value (i.e., loge(O) = 7.15) 
for T cff = 5650 K. We now point out that among our "best choice" sets of stellar parameters 
only one has a T e ff-value (T cS = 5450 K) which considerably deviates from that with 5650 
K. Such relatively low value for the effective temperature (i.e., T eS = 5450 K) causes a large 
disagreement between the NLTE Ol and [Ol]-based abundances amounting to ~ 0.34 dex. 
Moreover, as we have indicated above the set of parameters T e ff = 5450 K and [Fe/H]= —2.1 
produces an unpleasant scatter of the Fe il-based abundances around the mean value. In view 
of these factors we think reasonable to eliminate such particular set of stellar parameters. 

The remaining three sets of stellar parameters, i.e. T c g = 5670 K & [Fe/H]= —2.3, 
T eff = 5700 K & [Fe/H]= -2.3, and T cff = 5600 K & [Fe/H]= -2.0, correspond to the case 
of 1D+LTE modeling for both Fe I and Fe II lines, to 3D+LTE modeling of Fe II lines, and 
to 3D+NLTE modeling of Fe I lines, respectively. It is important to note that this range of 
effective temperatures is similar to the typical uncertainty in the determination of T e ff via the 
application of the infrared flux method. The small differences in the effective temperature 
implies a rather small oxygen abundance variation. In particular, our 1D+NLTE oxygen 
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abundance determination from the Ol IR lines using the former set gives loge(O) = 7.13, 
while the latter set leads to loge(O) = 7.20. For the [Ol] line we get loge(O) = 7.16 and 
loge(O) = 7.11, respectively. 

We must consider now how the dependencies shown for the oxygen lines in the rhs 
panels of Fig. 5 are going to be modified due to the effect of granulation inhomogeneities. 
As mentioned above, a major obstacle to investigate this point is the lack of a suitable grid 
of 3D models. In any case, this 3D issue concerns especially the [O i] line given its significant 
sensitivity to granulation inhomogeneities. According to Nissen et al. (2002) the main effect 
is due to differences in the continuum opacities between the ID and 3D models. Since the 
opacity is provided primarily by H~ atoms the strength of the weak [O i] line turns out to 
be inversely proportional to electron density. In the 3D hydrodynamical model of the star 
HD140283, which has been obtained assuming T cff = 5700 K and [Fe/H]= -2.5, the deficit of 
electrons is particularly large in the cool surface layers of the granular regions, while it is less 
pronounced in both the intergranular plasma and in the ID model. As a result, the spatially 
averaged equivalent width of the [O i] line in such a 3D model increases and, consequently, the 
derived oxygen abundance decreases in comparison with the ID modeling result. Obviously, 
the 3D correction obtained in Sect. 6.2 for this forbidden line assuming that the stellar 
parameters are T e g- = 5700 K and [Fe/H]= —2.5 cannot be used for the 3D+NLTE case with 
lower effective temperature (T cS = 5600 K) and higher metalicity ([Fe/H]= —2.0). 

As pointed out by Asplund et al. (1999) the deficit in the number density of electrons 
depends mainly on the metalicity. 3D hydrodynamical simulations for higher [Fe/H] values 
will have larger surface temperatures and electron densities. Thus, one can expect that at 
[Fe/H]= —2 the 3D abundance correction for the [Ol] line has to be smaller. Obviously, the 
exact value of such corrections will have to await the development of a new 3D model for 
the revised stellar parameters. It should be noted that such a sensitivity to 3D effects may 
be assumed also for the UV lines of OH. However, the impact on 3D abundance corrections 
when changing the metallicity from —2.5 to —2.0 seems to be rather small if we suppose 
that the results obtained by Asplund & Garcia Perez (2001) for dwarfs can be applied to the 
subgiant HD 140283 (see their Table 2, which shows that the 3D effects on the UV lines of 
OH is ~ —0.5 at [Fe/H]= —2 and ~ —0.6 at [Fe/H]= —3). The additional conclusion that 
follows from their Table 2 is a high sensitivity of the oxygen abundances derived from the 
OH molecules to the metallicity of the star. As a result, a 3D model with higher metallicity 
([Fe/H]= —2) might bring the oxygen abundances derived from molecules (loge(O) = 7.20 
for T cff 5800) to closer agreement with the abundance derived from the Ol triplet ( 
loge(O) « 7.20 for T cS = 5600 K.) 

An extra problem is related with observations of the [Ol] line in the star HD140283. 
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According to Nissen et al. (2002) this weak line (W ~ 0.5 mA) is blended with a stronger H 2 
line (W ~ 1 mA). These authors estimated that the available equivalent width measurements 
have an uncertainty of ±0.2 — ±0.3 mA . Such an uncertainty is very large, because it is 
about 50% of the observed equivalent width, and implies an uncertainty of ±0.13 — ±0.19 dex 
in the derived oxygen abundance (Fig. 5, middle right panel). 

Alternatively, in a first approximation we could neglect the sensitivity of the O I lines 
to the granulation inhomogeneities. Furthermore, the observation errors for these lines in 
HD 140283 are relatively small. For example, the observed equivalent widths of the Ol 
triplet that have been published in the literature (Abia & Rebolo 1989; Tomkin et al. 1992; 
Nissen et al. 2002; G. Israelian, private communication), and which we have used in our 
investigation, differ by about ±1 mA. Their impact on the derived oxygen abundance lies in 
the range ±0.06 dex (Fig. 5, middle right panel). It is clear that the O I lines seem to be the 
most reliable lines for deriving the oxygen abundance in the star HD140283. 

Therefore, when the stellar parameters are revised neglecting 3D and NLTE effects for 
the iron lines, and when the oxygen abundance derived from the O I lines is corrected for 
NLTE effects, then the mean value of the oxygen-to-iron ratio seems to be [0/Fe]~ 0.7 
at [Fe/H] =— 2.3. This value is in good agreement with the quasi-linear trend of [O/Fe] vs. 
[Fe/H] obtained by Nissen et al. (2002) using standard ID modeling. However, the key point 
to keep in mind is that a 3D±NLTE modeling of the O I IR triplet (using the 3D±NLTE 
stellar parameters derived from the Fel lines) would give [0/Fe]~ 0.5 at [Fe/H] =— 2.0. 
Obviously, the analysis of a single metal-poor star is not enough to opt for one of the 
following two possibilities: (a) a linear rise in the [O/Fe] ratio vs. metallicity (see, e.g., Fig. 
11 of Israelian et al. 2001) or (b) a plateau between 0.4 and 0.5 dex in the metallicity range 
-3<[Fe/H]<-l. 3 However, we think that our investigation indicates the "road to be taken" 
towards the resolution of the problem. 



3 We point out that with the stellar parameters used by Asplund et al. (1999) for the 3D modelling of the 
star HD140283 (i.e., T off = 5700 and [Fe/H]=-2.5) the NLTE oxygen abundance (loge(O) = 7.08 ± 0.06) 
derived from the Ol IR triplet and the NLTE Fei-based iron abundance (loge(Fe) = 5.77 ± 0.186) leads to 
a small oxygen-to-iron ratio ([0/Fe]~ 0.1 dex). Moreover, there is no overabundance of oxygen at all if one 
uses the [O i]-based abundance. It is also of interest to mention that in our preliminary study (see Shchukina 
et al. 2003) we obtained a slightly larger value for [O/Fe] because we used a smaller number of Fel lines. 
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9. Conclusions 

We have derived the solar oxygen abundance from the O I IR triplet using Asplund's 
et al. (2000) 3D hydrodynamical model of the solar photosphere and a NLTE modeling 
approach. We find loge(O) = 8.70±0.06, which is in good agreement with the determinations 
of Allende Prieto et al. (2001), Shchukina et al. (2003) and Asplund et al. (2004). Actually, 
if for the O I IR triplet we use the same oscillator strengths adopted by Asplund et al. (2004) 
we then obtain a value very similar to that reported by them (i.e., loge(O) = 8.65). 

In agreement with previous investigations, we confirm that LTE synthesis of Fe I lines 
in the 3D hydrodynamical model of this star gives flux profiles significantly different from 
those obtained in a ID model for the same stellar parameters. However, when the LTE 
approximation is used for the synthesis of Fe II lines we find negligible differences between 
the ID and 3D cases. 

The differences between the iron abundances obtained from Fel lines when assuming 
LTE in the ID and 3D models of the star HD140283 turn out to depend on the lower exci- 
tation potential of the spectral line under consideration. The NLTE abundance corrections 
for Fel lines are significant, especially for the low-excitation lines in the 3D model. 

We have three particularly important results concerning the determination of the iron 
abundance in metal poor stars like HD 140283: 

(1) If NLTE effects are taken into account when synthetizing the Fel lines in the ID 
and 3D models of this star, then the derived iron abundance turns out to be practically the 
same in both cases. 

(2) Contrary to a generally accepted belief, the Fell lines turn out to be significantly 
affected by NLTE effects. In particular, a full NLTE modeling should be carried out for the 
Fe II lines with upper levels of odd parity starting from z 6 D°, and especially the weaker lines. 

(3) The iron abundance of the star HD 140283 is estimated to be ~ 0.5 dex higher than 
previously thought. 

For the O I IR triplet the NLTE correction of the oxygen abundance in the ID and 3D 
models is practically the same and less than 0.2 dex. Such a NLTE correction is negligible 
for the [O i] line. We have confirmed that the mean value of the oxygen abundance derived 
from the Ol IR lines is insensitive to 3D effects, while such effects reduce the abundance 
derived from the [Ol] line by ~ 0.2 dex. 

With the here presented NLTE calculations for Fe and O, both the ID and 3D models 
lead to inconsistencies in the iron and oxygen abundance determinations. The discrepancies 
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in the derived iron and oxygen abundances cannot be removed by taking into account NLTE 
and/or 3D effects. This result has led us to investigate whether these discrepancies can be 
resolved by a modification of the stellar parameters of the star HD 140283. We have shown 
that previous studies may have underestimated the metallicity of this star and overestimated 
its effective temperature. We find [Fe/H]~ —2.0 (instead of [Fe/H]~ —2.5) and T e g = 5600 
(instead of T c g = 5700). With these new stellar parameters the iron and oxygen abundances 
in the star HD140283 would be loge(Fe) = 5.54±0.14 and loge(O) = 7.20±0.06, respectively. 
Taking into account the low value for the solar oxygen abundance (i.e., loge (O)=8.7O±O.O6) 
and our previously determined value of the solar iron abundance (i.e., loge Q (Fe)=7.50±0.10; 
see Shchukina & Trujillo Bueno 2001) we find that the oxygen-to-iron abundance ratio is 
[O/Fe] « 0.5 at [Fe/H]=-2. 

Besides the reported NLTE effects for the Fe II lines, our main conclusion here is that the 
metalicity of this type of metal-poor stars might well be significantly larger than previously 
thought. Obviously, the analysis of a single metal-poor star is not enough to fully resolve 
the puzzling behavior of the O/Fe ratio in metal-poor stars. However, we think that the 
present investigation at least indicates the road to be taken. 

We are very grateful to Ramon Garcia Lopez and Garik Israelian for allowing us to use 
their stellar observations of iron and oxygen lines. We are also grateful to I. Vasiljeva for 
helping us with the radiative transfer computations during the early stages of this project. 
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Fig. 1. — Departure coefficients {(3) vs. height in the 3D atmospheric model of the metal-poor 
star HD140283. Left panel: Representative "granular" point. Right panel: Representative 
"intergranular" point. The solid lines indicate the /^-coefficients for the Fel levels with exci- 
tation potentials below 5 eV. The dashed and dashed-dotted lines refer to the /9-coefficients 
of the lower and upper levels of the Fell lines in Table 1, respectively. The thick lines cor- 
respond to the Fel ground level a 5 D, to the Fell ground level a 6 D, and to the Fell excited 
level of odd parity z 6 D°, respectively. Arrows mark the representative atmospheric heights 
below which one of the Fel lines (A5429.699A) and the Fell lines from Table 1 originate in 
the NLTE and LTE cases. Note that such NLTE and LTE heights are nearly the same for 
the Fell lines. 
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Fig. 2. — The emergent disk-center intensity profiles of the Fel line A5429.699A calculated 
in the 3D hydrodynamical model of the metal-poor star HD140283. Left panels correspond 
to the NLTE case, while the righ panels to LTE. Top panels: the wavelength variation of the 
emergent intensity that one would be able to extract at each point along the spectrograph's 
slit if the 'surface' of the star HD140283 could be spatially resolved. The "wiggles" are caused 
by the Doppler shifts associated with the velocity flows of the stellar granulation pattern. 
Bright strips correspond to granules, while dark ones to intergranules. The hypothetical slit 
is located along the position Y=24 in the 3D model. Intensities are given in absolute energy 
units (erg cm _2 s _1 ster" 1 Hz~ 1 ). Black lines show variations of the line equivalent widths W 
along the slit. Bottom panels: The spatially resolved emergent intensity profiles at each 
of the 50 surface gridpoints considered. The thick solid lines with filled circles show the 
resulting spatially averaged intensity. 
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Fig. 3. — The flux profiles of the weak Fell line at 6247.56 A (top panels) and of the mod- 
erately strong Fell line at 5316.62 A (bottom panels) calculated in the 3D hydrodynamical 
model of the metal-poor star HD140283 (logg = 3.7, T cS = 5700 K and [Fe/H]= -2.5). The 
profiles were synthesized using the iron abundance corresponding to the metallicity loge=5.0. 
Left panels show the results for the NLTE case while the right panels refer to the LTE cal- 
culations. Individual thin curves show the computed emergent fluxes for the "granular" and 
"intergranular" models corresponding to the 50 spatial grid points situated along one of the 
F-directions in the 3D model atmosphere of the metal-poor subgiant HD140283. The thick 
solid line with filled circles in each of the panels shows the resulting spatially averaged flux. 
Fluxes are given in absolute energy units (erg cm _2 s _1 Hz _1 ). The profiles associated with 
the dark intergranular lanes of the stellar granulation have a redshift and a lower continuum 
flux. Bright granules are characterized by a higher continuum flux and a blue line shift. 
Note that the flux profile of the Fell line at 6247.56 A is observed in absorption (W > 0) 
while the synthesized NLTE flux profile turns out to be in emission (calculated W < 0). 
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Fig. 4. — The iron abundances derived from Fel lines (filled circles) and from Fell lines 
(stars) vs. the lower level excitation potential of the chosen spectral lines. Left panels 
show the results for the 3D model atmosphere of the star HD 140283 while the right panels 
for the ID model. Top panels: NLTE. Bottom panels: LTE. The four weak Fe II lines 
(AA6149.238, 6247.562, 6432.654, 6456.376 mA) used by Nissen et al. (2002) are marked with 
squares. The dashed and dotted horizontal lines indicate the mean abundances found from 
the set of Fe I and Fe II lines, respectively. The corresponding values of the mean abundances 
and their standard deviation are given inside the panels. 
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Fig. 5. — The computed abundances vs. the metallicity (top), vs. the effective tempera- 
ture (middle), and vs. the stellar surface gravity (bottom) for the set of ID stellar model 
atmospheres. Left panels: results for iron, with the solid and dashed lines indicating the 
mean Fe I- and Fe n-based abundances, respectively. Right panels: results for oxygen, with 
the solid and dashed lines indicating the mean oxygen abundances derived from the infrared 
triplet A7772-5 A and the oxygen abundance obtained from the forbidden [O i] A6300 A line. 
Thick lines in each of the panels correspond to NLTE while thin lines to LTE. The standard 
deviations (cr) for the observed equivalent widths of the O I and [O i] lines are shown in the 
middle panel. It is also of interest to note that our LTE+1D calculations show that the lower 
the effective temperature the larger the discrepancy in the iron abundance derived from Fe 
I and Fe II lines (see the thin lines in the middle panel to the Ihs of the figure). 
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Table 1. Ol, Fel and Fell line list and observed equivalent widths (W) from the 
metal-poor star HD140283. The oscillator strengths for the O I IR triplet are from Kurucz 
& Bell (1995), while that for the forbidden [Ol] line is from Lambert (1978). The 
wavelengths and excitation energies of the Fe lines are from Moore (1959). The Fe 
oscillator strengths were taken from Fuhr et al. (1988), except for Fel A4733.596 A and 
Fell AA5197.569, 5234.620, 6149.238, 6456.376 A which stem from Gurtovenko & Kostik 
(1989). The observed equivalent widths of all the Fe I lines and of 11 Fe II lines have been 
kindly provided to us by R. Garcia Lopez (private communication), while the equivalent 
widths of the remaining Fe II lines have been taken from Table 3 of Nissen et al. (2002). 
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Table 1 — Continued 
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Table 2. The solar oxygen and iron abundances obtained by different authors using 
Asplund's et al. (2000) 3D hydrodynamical model atmosphere and the ID photospheric 
model of Holweger & Miiller (1974). The given uncertainties reflect only the line-to-line 
scatter. The results are based on permitted O I IR lines, forbidden [O i] lines, OH 
vibration-rotation lines, OH pure rotation lines and lines of Fe I. 
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Table 3. The oxygen and iron abundances for HD140283 obtained by different authors 
using Asplund's et al. (1999) 3D hydrodynamical model atmosphere and ID hydrostatic 
models. The given uncertainties reflect the line-to-line scatter. In general, the abundances 
are based on different selection of lines, ^/-values and/or equivalent widths. The results 
given in the table are based on permitted O I IR lines, forbidden [O i] lines, Fe I and Fe II 

lines. 



Lines Approach Present study Nissen et al. (2002) Asplund et al. (1999) 
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Table 4. Stellar parameters of the star HD140283 
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[Fe/H] 
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